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Magnetic quasi-particles in hot QCD predict Casimir scaling for the 

tension a of spatial Wilson loops. Recently this scaling has been verified 

Z^ , on the lattice to 1-2% precision. We review how spatial loops tell us about 

the nature of the plasma state. The strong first order nature for more 
V^ I than five colours and the fast drop of the topological susceptibilty at Tc 

^D ■ found in simulations of hot quenched QCD are consistent with a quasi- 

"nI ' particle picture valid down to the critical temperature. The magnetic quasi- 

particles will Bose condense in the ground state below Tc, and may behave 
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r"| ' like a dilute Bose gas down to T=0. The model fixes the diluteness as 

CIh! "'/'^m- This is a small number for all temperatures. 
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1. Introduction 



Plasmas are ubiquitous in the universe. Most of visible matter is in this 
C^ I state, in which atoms are fully ionized. A plasma we know quite well is 

the sun. It sustains the thermonuclear power generation, that one tries to 
imitate in Tokamaks. Large static magnetic fields show up on the sun's 
surface as sunspots. They tell us that there is no magnetic screening effect. 
Of course static long range electric fields are absent, due to Debye screening. 

Today we are probing into a new plasma state, that of QCD, at RHIC. 
This plasma state features again Debye screening in the electric sector. But 
it is quite different in the magnetic sector. It has magnetic screening for 
all temperatures. At low temperatures 't Hooft [1] argued for a magnetic 
screening length. Such screening was expected on the basis of the dual super- 
conductor analogy for QCD. At high temperatures the need for a magnetic 
screening length was recognized also in the early days [2] of hot QCD. 



Based in part on lectures at the XLII Cracow School of Theoretical Physics, Zakopane, 
Poland May 30- June 8, 2003. The present version covers more material, and corrects 
some errors in the version published in the Proceedings of the School. 

(1) 
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It was based on the simple observation, that the static magnetic sector of 
QCD is three-dimensional Yang-Mills theory and that this theory contains 
a mass gap. 

This mass gap is, alas, not computable in terms of the small gauge cou- 
pling constant. But it is accessible by lattice simulations. On the basis of 
these simulations we find that the ratio of the square of the magnetic screen- 
ing mass and of the string tension of the square of the magnetic screening 
mass is a small number, about 1/20 with small corrections of 0{1/N'^). 

Some time ago a model for the three dimensional Yang-Mills theory 
was proposed [4], in which precisely this ratio was supposed to be small. 
The model is based on the hypothesis that the partition function can be 
approximately computed by thinking of the theory as Ibeing a dilute 3D 
gas of lumps with size equal to the magnetic screening length. The lumps 
are non-perturbative quantities in terms of the 3D gluons. But they are 
supposed to be in an adjoint magnetic SU(A^) multiplet. 

The merit of this model lies in simple consequences for the tension of 
Wilson loops with N-allity k. Their k dependence comes in through a simple 
multiplicity factor k{N — k). Since then this dependence has been verified to 
1-2% by lattice simulations. Apart from the multiplicity factor the model 
predicts the tension a of the loop to be simply the product of the magnetic 
screening /m and the density nyi of the lumps: 



<7~^M"-M- (1-1) 



Multiplying this relation with /^ gives the desired diluteness /^j^m- This 
diluteness equals l^j^cj and as remarked above, we know from lattice simula- 
tions this is a small number. Thus we have an a posteriori justification for 
the model. 

In the next section |21 we recall some basic facts and numerology for the 
Debye screening of electric charges, and how to simulate them on the lattice. 
We introduce the electric flux loop in a pedagogical manner. 

In section |31 we formulate the effective actions at high temperature. 
Those readers mostly interested in the quasi-particle picture can skip this 
section and the next on the critical behaviour and universality. In the 
next section 121 we pursue the same but now for the magnetic screening. 

Finally we pass to the subject of magnetic quasi-particles in section [HI 

The last section contains a brief discussion of the first order nature of the 
transtion at large number of colours, conclusions and prospects. Appendices 
do elaborate on some intricacies, especially on the use of a simple Stokes 
theorem. 
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2. The Debye mass 

Let us consider some gas heated at a temperature well above the ioniza- 
tion energy. The electrons and ions are then moving more or less indepen- 
dently from each other. For simplicity we take the ions to have the opposite 
charge of the electrons. 

Since the ions are much heavier than the electrons one can consider 
them to be a charged background with a density \e\n. If you immerse a 
heavy point charge Q in this medium its Coulomb interaction changes the 
density n{x) of the electrons around it. The ions will remain unperturbed. 
As all of you know this gives rise to screening of the Coulomb law. The 
argument is purely classical. The Poisson equation in the presence of the 
ionized electrons with charge e reads: 

A^o = en{x) - ern - Q5{x) . (2.1) 

The variation in the electron density is due to the variation in the en- 
ergy cAq of electrons in the field Aq , and hence in the corresponding Boltz- 
mann factor exp {eAQ/T){ setting Boltzmann's constant equal to 1). Putting 
Aq{x = oo) = we get n(x) = nexp {eAQ{x)/T). The Poisson equation be- 
comes in the linear approximation in the energy: 

2 
P 77 

AAo = —Ao-QS{x). (2.2) 

Eq. (|2.2j) is solved by a Yukawa potential, ^o = Q '^^'^4~™°^ , with the 
Debye screening mass mn: 

ml = —. (2.3) 

The screening length /d is the inverse of the screening mass rriD- Its 
raison d'etre is statistical, due to the Boltzmann factor. So we expect the 
screening to involve many electrons. And this is precisely what Eq. 1)2. 3(1 
tells us: in a sphere of radius /d we have 

lln= r3/V(e3ni/2) (2.4) 

electrons. 

To estimate this number we introduce the atomic radius r^ to split up 
the expression in two dimensionless factors: 

lln={T/{e^/ra)fl\rlnr'l^ (2.5) 

The first dimensionless factor is large because T is larger than the ioniza- 
tion energy e^ jr^- The second dimensionless factor is 0(1) or larger because 
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the number of electrons inside the atomic radius is certainly not larger in 
the ionized state. 

Typically, for the sun's corona T/e^ is about 10^^^ the atomic scale and 
the number of electrons in the Debye sphere is 10^. This condition is called 
the statistical screening or plasma condition. 

It is useful to note that the thermal de Broglie wavelength A^ is small 
with respect to the interparticle distance. In fact the thermal wave length is 
on the order of the atomic radius and hence much smaller than the interpar- 
ticle distance. The thermal wave length, just above ionization temperature, 
but well below pair creation temperatures is, in units of the atomic radius : 

* rp (-1/2) 

AbK = /i/(V27rmTr„ = V2^ ^ - < 0(1) (2.6) 

We used that w, the electron velocity around the nucleus, is a times the 
velocity of light, and that the Bohr condition fi/ravTa is on the order of 1. 
This estimate is valid for temperatures above ionisation but well be- 
low pair-creation scales. So the physics of the electron plasma just above 
ionization is classical. This is no longer true when pair creation becomes 
important. 

2.1. A Gedanken experiment 

It is amusing to do the following Gedanken experiment. It will provide 
us with a characterisation of the ionised state, or order parameter. Suppose 
we want to compute the electric flux going through some large (with respect 
to the atomic size) closed loop L with area A(L). Normalize the flux <l> = 
JgdS ■ E by the electron charge e and define : 

ViL) = expi2TT^/e. (2.7) 

Of course, at T below the ionization temperature no flux would be de- 
tected by the loop, because there are only neutral atoms moving through 
the loop. 

Let us now raise the temperature above Tjonisation- What will happen? 
Both electrons and ions are screened. For simplicity we will take the ions 
to have the opposite of one electron charge. 

We are going to make the following simplification. The charged particles 
are supposed to shine their flux through the loop if they are within distance 
If) from the minimal area of the loop. This defines a slab of thickness 21 jy. 

Of course if we plot \'P/e\ as function of the distance of the particle to the 
loop you find an exponential curve with the maximum 1/2 at zero distance. 
For the sake of the argument we will replace that curve by a theta function 
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of height 1/2 and width 2/d- If one wants to do better one has to deal with 
infinitesinially thin slabs, and integrate over the thickness, and we will do 
so in Appendix A. The result is parametrically the same as the one we will 
derive keeping the simple minded method. 

Then one electron (ion) on the down side of the loop will contribute 
+l/2(— 1/2) to the flux, and with opposite sign if on the up side of the 
loop. That is: l^(-^)|one charge = " 1- This result is independent of the sign 
of the charge! The plasma is overall neutral, the loop is sensitive only to 
charge fluctuations. For / charges inside the slab the flux adds linearly and 
we find: 

Vim = (-!)'• (2.8) 

Assuming that all charges move independently, the average of the flux 
loop V{L) is determined by the probability P{1) that / electrons (ions) 
are present in the slab of thickness 2/d around the area spanned by the 
loop. Taking for P{1) the Poisson distribution jt(0^ ^^P ~^~ ^i^ the average 
number of electrons (ions) in the slab- we find for the thermal average of 
the loop: 

{V{L))t = J2mV{L)\i = J2Pm-y = exp-4[. (2.9) 

I I 

Now I = A(L)2l-£,n(T) , so the electric flux loop obeys an area law 
exp —p(T)A{L), with a tension 

p{T) = 8hniT). (2.10) 

Note that the tension in units of the Debye mass, l^p, equals the number of 
electrons in a Debye volume 8/£,n(r), because of ea. ()2.1U() . We saw already 
that this number was huge. 

A caveat behooves the word "tension" . On the minimal area the density 
of electrons is not changing, only near the perimeter. Nor is the free energy 
density. So the word tension is not the correct one, though we will continue 
to use it. 

Another aspect of this tension is that no quantum theory came into the 
calculation. The only ingredient that takes into account the dynamics of the 
plasma is the screening length Id and the statistical independence through 
use of the Poisson distribution. 

In conclusion the behaviour of the loop is very different for the ionised 
state, as compared to the de-ionised state. It behaves with an area law 
above ionization temperature. Un-ionized atoms can at most cause perime- 
ter effects. 
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2.2. A closer look at the flux loop 

Untill now only classical physics went into the calculation of the flux 
loop average. 

If we assume for the electron density n the Fermi-Dirac distribution our 
calculation is semi-classical. But only at a temperature where relativistic 
effects are starting to operate (T ~ 2me) the difference between the Poisson 
distribution and Fermi-Dirac becomes important. For temperatures, where 
electrons and positrons become easily pair produced we need a fully-hedged 
field theory. As we will see at the ned of this section we can compute p in a 
WKB approximation. Such approximations are of 0{1/ g) and are therefore 
large. 

And our Gedanken experiment takes in field theory a perhaps surprising 
form. The flux /^ dS ■ E becomes now an operator, E being the canonical 
conjugate of the vector potential A, with [£^fc,yl;] = jhS^^i, the Kronecker 
delta standing for a delta function of the space arguments as well. 

The loop acts on the state space, and in particular we are interested in 
what it does to physical states. What are physical states? First we should 
observe that at the ambient temperature T > 2me the nuclei start to come 
apart. So we have a field theory of electrons, positrons, and nucleons. A 
physical state is one where a given state contaning a set of particles is made 
gauge invariant. How to do this can be seen from a gauge transformation 
on the photon field: 

Ak -^ Au = Ak -\ — dkUJ = Ak exp iujdk exp —ico (2-11) 

e le 

and for the matter fields: 

i)^e^p{ioj)tp. (2.12) 

So the gauge group is U{1), and lo is the angle characterizing the U(l) 
phase expiu lying on the unit circle. Now a physical state is produced 
by averaging a given state over all possible gauge configurations uj. "All" 
excludes those configurations where you can find a closed curve along which 
CO winds around the unit circle once or more. That is to say, all admitted 
gauge transformations constraining the physical Hilbert space are in the 
trivial homotopy class Hi of the unit circle U{1) ^. 

To understand the action of the flux operator V{L) we consider a gauge 
transformation uJi{L). The suffix indicates a discontinuity of 2tt when we 
go with the gauge transformation through the area spanned by the loop L. 



^ Had we constrained the physical Hilbertspace even more, by including transforma- 
tions in the non-trivial homotopy classes, the loop V{L) would have had no physical 
effect. 
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One can take for the gauge function half the solid angle spanned by the 
loop. This makes a jump of 2tt at the surface, and vanishes at infinity. 
Gauss' theorem tells us then: 

f dS.E = ^ / dVV{EuJi{L)) (2.13) 

with the volume integral over all space. This formula is true for any ui{L) 
with discontinuity 2tt on the surface S{L), and which falls off sufficiently fast 
at spatial infinity, so that there we are allowed to drop the surface integral. 
The r.h.s. can be rewritten with the help of the charge density operator ej'o 
as: 



■ J dVV{Eui{L)) = ^J dV[{V.E - eJo)ui{L) + (^.V + eJo)u;i{L)] 



2tt 

(2.14) 

So the flux operator can be written on the physical subspace ( where 
V.-B — ejo = 0) as a gauge transformation with the singular gauge func- 
tion. 



y(L)=exp{i- f dV{E.V + ejo)uji}. (2.15) 

According to this formula V{L) creates a magnetic flux line along the 
closed loop L. Everywhere else it is just a regular gauge transformation. So 
its physical effect is localized on the loop L. Hence the volume integral in 
eq. H2.15() reduces to a line integral along L, that we can write formally as: 

V{L) = expi- I dV{E.V + ejo)wi = expi f dl.fl. (2.16) 

The flux dl.fl has strength / ds.jVuJi = -^, where 7 is any closed curve 
encircling L once. This means that a Wilsonloop: 

W{-i) = expi| i ds.A (2.17) 

will commute with the electric flux loop. And this is borne out by using the 
canonical commutation relations. 

It is amusing that in the process we have uncovered a dual Stokes the- 
orem that relates V{L) as line integral of this magnetic flux to V{L) as 
surface integral of the electric field strength. It is only valid in the physical 
subspace. 

Loops with an integer multiple k of unit flux are possible too. They 
will be denoted by Vfc(L) and correspond to the n-th homotopy class of the 
homotopy group Iii{U{l)). Our "dual Stokes theorem" then tells us that 
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27r /■ - - 
Vfe(L) = expi— A; / dS.E. (2.18) 



This concludes our simple-minded exploration of the flux properties of 
the plasma. The question is: can we do better? To answer that question 
we have to use a more sophisticated version of our screened electron gas. 

2.3. The Coulomb gas 

One reproach one can make to the model in the previous section is its 
treatment of the thermal fluctuations. The only way they are built in is 
through the Debye mass and the Poisson distribution. To take them more 
fully into account we are going to analyse the three dimensional Coulomb 
gas. This model still ignores the time dependent fluctuations in QED. 

The Coulomb gas has as degrees of freedom the density n{x) = J2£- ni6{x- 
Xi). Ui is the sum of all the charges in the point Xj. The interaction is by 
means of the Coulomb propagator: 

D{x-y) = {x\^\y) (2.19) 

The partition function is then given by: 

Z = 5^ exp -|- / ^n{x)D{x - y)n{y) - /i^ [n{xf. (2.20) 



n(x) 



The parameter fi has dimension of mass, and serves to control the charge 
density. 

We introduce a continuous parameter A{x) that decouples the densities 
in the partition function. In what follows we are going to simplify notation 
by replacing integrals by dots and dropping the position vectors. So we 
have: 

Z= f DAj2exp{--{{^)VA.VA)-in.A- n^n'^}. (2.21) 

Whatever the value of /i, after summing over n we find an effective 
potential V{A), periodic mod 2it in A: 

Z = J DAexp- J^{^(VAf + V{A). (2.22) 

Apart from periodic, the potential is even in A. So it can be written as 
a sum of cosines. For small /i only the cos A term appears. So for small A 
we get a mass term e'^fi^/T, so the A field is the screened scalar potential. 
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Fig. 1. The profile of the field A that minimizes the effective action in ea. ()2.22|l . 

Let us finally compute the average of the electric flux loop in this theory. 
Remember the flux loop can be formulated as any gauge transform in three 
dimensional space, where the gauge parameter to makes a jump of 2tt at the 
surface. The gauge parameter is nothing but our variable A(x)\ 

Thus, the average of the loop is determined from the effective action 
in ea. (|2.22|) . where A now depends only on the coordinate normal to the 
loop surface, vanishes at infinity and makes a jump of 27r at the surface-see 



<Vi{L) >= J DAexp-Ll J dz{^{d,Af + V{A)). (2.23) 

As the transverse size of the loop, Ltr — > oo we only need the minimum of 
the effective action. 

Since the potential V is periodic mod 27r, we can lift the branch on 
the right hand side of the loop area in fig.(l) over 2tt. Then we have a 
continuous profile starting at A = at — oo and ending at A = 2tt at +oo. 

This is the solution of the classical problem where the profile A(z) rolls 
from one maximum to the next one in V. Again we found an area law for 
the loop, the tension being 



pi = miuAiJ dz{^{d,Af + V{A))} 



(2.24) 
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Note that the potential V(A) does not depend on the coupUng e. One 
has to absorb the coupUng e into the z-coordinate, in order to balance kinetic 
and potential energy: the profile that minimizes the action must be of the 
form A{ez). A scaling argument then leads to 

pi = c-T'^. (2.25) 

The precise knowledge of V does not matter for the parametric form. 
The numerical constant c reflects the ratio of T and //, and the form of the 
potential. The reader recognizes the same parametric form as in the result 
ea. ()2.1U() from the quasi-particle picture. 

2.4- Summing up 

Let us conclude his long section with a short summary. The quasi- 
particle picture gives an area law for the loop, just from the Poisson distri- 
bution. A more sophisticated version, the Coulomb gas, gives an area law 
as well. The new feature is its semi-classical nature: it is due to tunneling 
through an effective potential barrier. 

If we want to extend the calculation to field theory, with time depen- 
dent fiuctuations, we have to worry about Aq, the variable in our effective 
potential, since it is no longer gauge invariant. This is naturally mended by 
replacing the static ^o by the phase of the Polyakov loop or Wilson line: 

rl/T 

P{Ao) = V exp ie / Ao(r, x))dT = exp iA. (2.26) 

So in the full gauge theory we are supposed to compute the effective 
potential as a function of this phase. This can be done in a saddle point 
approximation around the phase A. Once done, the calculation becomes 
the minimization of that action. 

To one loop order the potential V equals the logarithm of the deter- 
minant of the fluctuations around A. This determinant has been studied 
extensively in the litterature [13] and the result for the tension is paramet- 
rically the same as what quasi-particles gave us. In terms oi q = i^;^ and the 
rescaled z-variable z' = eTz one finds for the effective action : 

2'k'^T'^ r 2 11 1 

Seff = —^ J dz'{{d.,qf + -(- -{q + -)\l -\q+ -\f)) (2.27) 

The potential y(g) = |(^ -(g-h^)^(l-|g-h^|)^) is the logarithm of the 
fermionic determinant in the background q normalized by the determinant 
without q. 
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The minimum of this action is easily found by completing the square of 
the kinetic term: 



The minimum is realized when the square cancels. This gives us the 
equation of motion for the profile q. The remaining term is the actual value 
for pi and equals 27r^T^ / dq^jV{q) through change of variables. 

Rewrite this as Cflon, using the lowest order result l^ = ^-3— and 

n = 0.87. ..{■^:^) , and you have the same parametric result as that from the 
quasi-particle picture. 

Only the numerical factor c/ is quite a bit larger than in ea. H2.1U|) . 
because the time-dependent fluctuations have now been taken into account. 



3. Effective field theories at high temperature 

With what we have learnt from the static approximation in mind we 
turn to QCD at high T. 

Any field theory in equilibrium at non-zero temperature can be formu- 
lated as a Euclidean path integral. The time direction in that integral is 
periodic mod 1/T for bosons, and anti-periodic for fermions. For the sta- 
tistical Gibbs sum one has 

TTp^ys ew-H/T = J DADqDqexp-^S{A,q,q), (3.1) 

where the gauge potentials A and the quark fields q are integrated over. 
Under gauge transformations we have: 

(3.2) 
(3.3) 

and covariant derivatives D{A) = d—ig[A, and V = d—igA. From the com- 
mutator of the partial derivatives follows F^y = d^Ay — dyA^ — ig[A^, Ay]. 
In the limit that T becomes small with respect to typical particle scales 
the time direction can be neglected. This is called dimensional reduction [6]. 
It can be formulated as a systematic approximation scheme using that QCD 
at high temperature T has a small running coupling g{T). The inverse 
propagator of a boson is proportional to {2imT)'^ +p^. Here n takes on 
integer values. For a fermion n is replaced by n -|- 1/2. As T becomes large 
one can integrate all heavy T-modes ( "hard" modes) and stay with a three 



A - 


-^ QAQ^ - 


-—ftdn 


q - 


-^ Qq 





12 C.P. KORTHALS AlTES 



dimensional theory in terms of only the static bosonic modes with n 
For QCD this 3d Lagrangian reads: 



/:e = TT{D{A)AQf + mlTiAl + \^{TT{Al)f 

1 
2 



+AE(Tr(^o)' - ^O^^Alf) + Tr4 + 5C^ . (3.4) 



This is called the electrostatic Lagrangian. 

The last term contains higher powers of Aq and of the covariant deriva- 
tive D{A) = d — igE[A. Neglecting it means one neglects 0{g'^) in the 
correlations you compute with the first six terms. The parameters in this 
Lagrangian are computed from the corresponding n-point functions in one 
and two loop accuracy in terms of g| = g'^{T)T. Higher loop order adds only 
to the accuracy if one takes into account (5£e- The electrostatic coupling 
g^ is to one loop order in terms of Ajjg (no flavours): 



glN _ 2Aw 



6.742.T\ ■ ^'^•^) 



T lllog(^i^ 

MS 

The subtraction was chosen to minimize the one loop effects [5]. 

We have swept one problem under the rug. When integrating the hard 
modes we have to admit a lower cut-off tIe, in between the scale T and 
the electrostatic scale gT. In principle the parameters will depend on this 
cut-off. 

One expects that we have the same picture as before: above the decon- 
fining temperature T^ ~ 200 MeV) we have a gas of "ions", the quarks, and 
of "electrons" , the gluons. There is screening as before, as witnessed by the 
mass parameter m^ = -^g^T^ in electrostatic Lagrangian, Eq. 1)3.4(1 . 

This constitutes the Stephan-Boltzmann picture of QCD and interac- 
tions between gluons and quarks describe deviations from this free quasi- 
particle system. 

Specific to QCD is that there is not only an electrostatic scale set by 
the Debye mass. We can integrate out in electrostatic Lagrangian all de- 
grees of freedom corresponding to the Debye scale. That leaves us with the 
magnetostatic Lagrangian: 

Cm = Tr4 + 6Cm (3.6) 

with a magnetostatic gauge coupling g^^ in Fij . 

Here an ultra-violet cut-off tIm is needed. It separates the electrostatic 
scale gT and the magnetic scale g^T. 

Terms with higher order covariant derivatives are contained in 6Cm- 
They are needed when we want an accuracy of 0{g^). 
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The magnetic action gives a non-perturbative theory. In calculating a 
Green's function with a typical external momentum p we find as dimen- 
sionless parameter g^ji/p which is 0(1) if the momentum is the magnetic 
scale. In particular calculation of the free energy in this theory will give 
for dimensional reasons ((7m)^ times a non-perturbative constant. That is 
the contribution one expects from a 4-loop diagram. As for the Green's 
function all higher loops are of the same order. 

Still we can compute a series in the small coupling g{T). Only the 
coefficients are non-perturbative. 

A useful comment on the parameters in the two effective actions is in 
order ^. Their determination is perturbative. Truncation of higher deriva- 
tive terms in the electro- (magneto-)static actions meant dropping terms of 
relative order 0{g'^) {0{g^)). This means we need to calculate the param- 
eters to these precisions. In particular we need gE to one loop order as in 
ea. (|3.5j) . and gu to two loop order in terms of gE and mg. The latter rela- 
tion between magnetic and electric scales by integrating the out the electric 
scales from the electrostatic action and matching with the magnetostatic 
action has been computed and gives [20] 

9lN = glNil - ^^ - ^/^? + 0{g')). (3.7) 

48 ■nniE 4dUo TTniE 

As niE is of order g the dimensionless expansion parameter is 0{g). 

For asymptotically large temperatures such a picture is indeed accurate. 
But asymptotic means temperatures about lO^Tc, well above the electro- 
weak scale, far away from where RHIC physics operates . It turns out, 
however, that the accuracy varies with the observable. 

In this subsection we will take the pressure as an example, how the 
expansion works. To put the calculation of the contributions of order higher 
than three in perspective and to see how the different scales come in, we 
recall once more the hierarchy of scales, cut-offs A and reduced actions 
needed to compute the pressure: 

T > tIe > ^T > TIm > c/^T . 

The pressure is normalized bypo = Pstefan-Boitzmann and consists of three 
parts: 

P 

— = Ph + PE + PM ■ 
PO 

The hard modes are cut-off in the infrared by tIe and equal ph- Schemati- 
cally we get: 

Pi, = l+g^ + g^ log --r + 9^ + 9^log-r- + g^ + ... . 

/IE ^E 



Brought up in conversations with M. Laine. 
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All powers of the coupling are even, since infrared divergencies are cut- 
off by A-E- The short distance scales (larger than T) are absorbed in the 
running coupling, Eq. 1)3. 5|) . The cut-off /Ie appears only in logarithms. 
The electric mode contributions are computed with Ce and give pE'- 

^E = 5'+/log^ + <7^+/ + /log^ + /log^ + / + .... 

Note the odd powers in g. They come in because the electric mass gT 
comes in through propagators from the electrostatic action, Eq. ()3.4|) . For 
example, the contribution from the scalar potential ^o gives the first term 
in pe: 

79^W^°S(A: +mE) = r 

V'T^r 167r2 

The dominant cubic term was computed in Eq. (|3.8() . We can expect 
logarithms of the two ratios of the three scales tue, Ae in the electrostatic 
action and Ay^. 

Finally the magnetic contribution is computed with Cm'- 



- ^(iV^ - 1) / 7^ log(P + ml) = ^ml . (3.8) 



6i ^M , 6 , 

9m 

We only put in the obvious dependence on the parameters in the elec- 
trostatic and magnetostatic actions. There are three comments: 

• All terms shown are perturbatively calculable, except the last one in 

PM- 

• All perturbatively calculable terms have been computed [16], except 
for the g^ terms. In particular the log's are known by now [9]. 

• All dependence on the cut-offs cancels, as expected. 

This is dramatically illustrated by Fig. [2 You see on the left the lattice 
data for the pressure in units of its Stephan-Boltzmann value ^g-T^ plotted 
together with the known low order (up to 0{g^)) perturbative results. For 
the asymptotically large temperatures mentioned the quasi-particle picture 
is indeed accurate as the figure shows. The right panel shows how the 
prediction can improve, when the known [9] logarithmic contribution to the 
0{g^) coefficient is included, together with a guess for the non-perturbative 
part of the coefficient. 

For any reasonable T, say below 2 GeV, the coupling obeys g1{T)/ 
T < 2.5. This is about 30 times bigger then e^, so we may already surmise 
that low order perturbation theory will be far from accurate. 
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Fig. 2. Left: perturbative results at various orders for pure SU(3) gauge theory, 
including 0{g^) for an optimal constant. Right: the dependence of the 0{g^) 
result on the (not yet computed) constant, which contains both perturbative and 
non-perturbative contributions. The 4d lattice results are from [7]. From Ref. [9]. 

3.1. The Dehye mass and electric flux loop in QED and QCD 

Now we discuss the Debye screening in the QCD plasma. 
Let us put a probe charge in the plasma, say a very heavy quark. In 
Fig. El the exchange of a single gluon is shown, together with its multiply 




Fig. 3. A single gluon exchanged between two static test charges. 

inserted self-energy. Once we compute the self energy 77oo of the gluon, 
shown in Fig. 0J the resumed propagator Dqq becomes: 



Doom 



p' + nooip) 



(3.9) 



It is easy to see that only the hard modes contribute to the self energy. 
Hence one finds the one loop result for ?tt-| in the electrostatic La- 
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grangian. 



exp -mEr . . 

-Dool?') ~ • (3.10) 



To two loop order one finds that already at that order non-perturbative 
effects contribute. 

Hence a definition independent of perturbation theory is called for, and 
a natural candidate is the correlator between two heavy test charges: its 
fall-off as a function of distance gives us the screening mass. 

The test charge put in the plasma changes the free energy. This change 
can be expressed in terms of an expectation value of the thermal Wilson 
line: 

Ai^V lDAJ^Tvr{Ao)ew-S{A) _^^^^ ^^ 
^^^--T = jDAe^p-SiA) = ^^("^°)^ ^'-''^ 

where the thermal Wilson line is given by: 

T 

P{Ao{x)) = ^TTrexpigldTAo{x,T). (3.12) 



If the test charge is in the fundamental representation then so is ^o- 
The path ordering is defined by dividing the interval [0, l/T] into a large 
number N^- of bits of length At = j^ : 

V{Ao) = lim U{t = 0, At)U{t = At, 2At) . . .U (r = ^ - At, ^) . 

(3.13) 
^From a formal and from a computational point of view the correlator 
{P{Aq)P{Aq)'^) has two advantages over the calculation presented above 
using the scalar potential: 

• The correlator is gauge invariant. 

• The correlator can be evaluated non-perturbatively, i.e. on the lattice. 
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Both are needed for an accurate determination of the screening length 
in QCD. Perturbation theory is not enough, despite the small coupling g{T) 
for asymptotically large temperature {i.e. well above the electro- weak scale). 
In the confined phase the correlator {P{AQ[r))P[AQ[0))^) obeys an area law 
exp— 0"^', in the deconfined phase the area law is replaced by the Yukawa 
potential controlled by the Debye mass. 

There is a further advantage: correlators of gauge invariant operators 
will excite the levels of a fictitious Hamiltonian describing Yang-Mills dy- 
namics in a space with one periodic mod 1/T direction and two other infinite 
directions. The time conjugate to this Hamiltonian is now the direction of 
the correlation. 

Conserved quantum numbers are then, apart from those from the two 
dimensional rotation group, the usual discrete parities, charge conjugation 
C, parity P (now in 2D), and a new quantum number, that changes Aq — > 
— ^0) called i?-parity. The magnetic quasi-particles will condense in the 
ground state below Tc. 

Note that the Debye mass defined this way should coincide in one loop 
order with what we found before: tti-e in the electrostatic Lagrangian. So it 
is associated with i?-parity —1. Clearly, this corresponds to the imaginary 
part of the Wilson line. The real part excites R = +, P = C = + states. 
To wit: if the correlator {PP) between two like charges were zero, then the 
difference between correlators of imaginary and real parts would be zero. 
That would mean, in turn, that the masses controlling their decays would be 
degenerate. This is not the case. Two like screened charges are compatible 
on a torus. In the confined phase (no screening) their correlation is indeed 
zero. 



4. Z(N) symmetry, universality, and the order of the transition 

Let us move from the very high T region to the critical region near 
Tc- There is a symmetry due to invariance of the Yang-Mills action under 
gauge transformations that are not periodic in Euclidean time, but only 
periodic modulo a center group element expifc^/jy. I^v is the N x N unit 
matrix. So with k integer, the determinant is one. What is not invariant is 
the periodicity of fields in representations with non-zero A^-allity, like quark 
fields. So we will discard them for the moment. 

Now the Wilson line, Eq. ()3.12|) . under such a transformation is mul- 
tiplied by the Z(A^) phase factor expik^. So the probability to find the 
system with a given value for the Wilson line: 

E{P) ~ / DA6iP - P{Aoj ) exp -S{A) (4.1) 
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has the same value in P as in expik^P, because the measure stays the 
same, so does the action, only the argument of the delta function will change. 
This is a most useful manifestation of Z(A^) symmetry. It is the analogue 
of the Z symmetry we met in the case of QED. 

Note that the Wilson line is a scalar quantity in every point x. 

So it bears resemblance to a Z(A^) spin variable z{x) defined on a three 
dimensional lattice. If we endow this Z(A^) spin system with a nearest 
neighbour Z(A^) invariant action, we have a system that has a transition 
point where the spin system changes from disordered into ordered behaviour. 

There is now the hypothesis [25] that the transition of this spin system 
and that of the Yang-Mills system are in the same universality class. This 
is interesting because it relates critical behaviour of a rather simple system 
to that of our Yang-Mills system. 

For N = 2 and 3 this spin system is unique in that one can write down 
only one action for this spin system per link: 

Sn=2,3 = /3{z + z*) . (4.2) 

Here z is a shorthand for the product of the two spin variables at the end 
points. 

Indeed the transition is second (first) order for A^ = 2(3), and many 
studies have found that critical behaviour is identical [10] . 

However for A'^ = 4 the spin system is not unique. Let us parametrize 
the action per link like: 

SN=i = l3{{z + z*) + xz'^). (4.3) 

In Fig. [21 [3] the phase diagram of this theory is plotted. Only positive 
couplings are of interest to us. 

The VEV of the z spin corresponds to the VEV of the Wilson line 
P{Aq) [12], the VEV of z"^ to that of P(ylo)^). The region in between the 
two second order transition lines corresponds to the subgroup Z(2) already 
broken, but not yet Z(4). That would imply two Debye masses (not very 
natural from the plasma point of view), one corresponding to P{Aq) and 
still zero in that region. Another corresponding to P{Aq)'^ and already non 
zero in that region. 

But Nature has decided differently: in the gauge system the transition 
is first order [10,7], and, from the phase diagram, that corresponds to both 
order parameters jumping at the same time. This is what has been con- 
firmed [11] in gauge theory within errors. 

Recent data [26] show that the first order transition becomes stronger 
with increasing N. This is consistent with the idea that quasi-particles gov- 
ern the behaviour of the plasma from very high T down to just above the 
critical temperature. 
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Fig. 5. Schematic phase diagram of three-dimensional Z(4) spin model from 
Eq. H4.3|l on a simple cubic lattice, taken from Ref. [3] , where it was extracted from 
series analysis and Monte Carlo data. Dashed and solid lines indicate first and 
second-order transitions respectively. Dotted lines indicate cases where the nature 
of the transition has not been unambiguously determined. The phases are labeled 
disordered ((z) — (z^) = ); Baxter (ferromagnetic with (z), (z^) both non-zero); 



"(si)" (where (z ) is ferromagnetically ordered but 



0). 



4.1. Electric flux and the spatial 't Hooft loop 

The phenomenon of deconfinement involves the breaking of the electric 
flux tubes, and the appearance of quasi-particles, the gluons. This reflects 
itself in the change in the force law between test charges, discussed above. 

How does it manifest itself in other measurable quantities? A natural 
candidate is the spatial loop that measures the electric flux, as we discussed 
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in the first section. This loop is formed by a closed magnetic flux line, the 
't Hooft loop [14]. 

What will perspire [4, 18] is that the behaviour of this loop in the de- 
confined phase is again quite different from that in the confined phase. The 
quantitative behaviour of the loop at very high T can be computed along 
the same lines as in Section |21 This is what we will do below. 

We start with a definition of the loop as a magnetic flux loop, i.e. as 
a gauge transformation e'xp {iui\^{x)YiS) with a discontinuity expifc^ when 
going around the loop. Here Y^ = diag jj{N — k, N — k, ...,N — k, —k, ... — k) 
with N — k entries —k and k entries N — k, so that it generates the center 
group element: 

exp {i27rYk) = exp -ik— = 4- (4-4) 

uJi,{x) is half the solid angle defined by the loop. 
In the Hilbert space this operator reads: 



/■ 1 ^ ^ Yu 

Vfc(L) = expi / dx-TtEix) ■ Duji^ix) — . (4.5) 

J Q 2 



9 
A representation which has the same effect in the physical Hilbert space 



is: 

Vfc(L)=expi— / dSTTE{x)Yk. (4.6) 

9 J 

S{L) 

Using the canonical commutation relations you can check that 14 (L) and 
Vk{L) have the same effect on physical states. That is, they multiply Wilson 
loops with the center group factor if the latter intersects with S{L). More 
precisely, if the Wilsonloop Wi{L) is in a representation with N-allity /, i.e. 
built from / mod N fundamental quark representations, then the 't Hooft 
commutation relation reads: 

Vk{L)Wi{L')Vk{L)^ = exp{m(L,LO— I^K^O- (4-7) 

The integer n{L,L') equals or 1, depending on how many times (mod 
2) L' intersects the surface spanned by L. The reader can find more detail 
in ref. [35]. 

Alternatively you can verbatim follow the proof of the equality given in 
section 13 So again, eq. (|4.6|) is the dual Stokes version of Eq. (|4.5() . 

Are the matrices Y^ unique? The answer is that they can be replaced 
by ^YkO,\ with Q a regular gauge transformation, without changing the 
effect of the flux operator in the physical subspace^. But that is all the 

^ This reflects the fact that the flux operator transforms hke an adjoint under a gauge 
transformation. 
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freedom one has. Other choices, hke generators of subgroups of SU{N), 
lead to gauge non-invariant flux operators. 

4-2. Quasi-particle calculation of electric loop average 

In the confined phase the particles are colour neutral so will at most 
contribute a perimeter law. At least, that is what one would infer from the 
analogy with the QED flux loop in the previous section. However, for QCD 
we have to distinguish glueballs and bound states of quarks. A glueball is 
understood to be a closed fundamental electric flux loop. So, once close 
enough to the loop it will send an fundamental flux through the surface, 
multiplying the k-loop with a centergroup factor Zi. But since it is closed, 
the flux has to come back through the surface giving a factor zi*. These 
phases interfere, and thus the total effect is as if there were no charge. Only 
at the perimeter there is a net effect. 

On the other hand a meson (aqq pair) close to the loop has an effect: it 
multiplies the loop with a center group factor Zi or its complex conjugate, 
depending on its orientation. Now the phases do not interfere, and a surface 
effect arises. 

Let us now turn to the deconfined phase and repeat the quasi-particle 
argument for the area law: 

{Vk{L))=exp-pk{T)A{L). (4.8) 

The gluons are in the adjoint representation, so their Y^ charge follows 
from the differences of the diagonal elements of Y^. That means there are 
2k{N — k) gluons with charge ±(7, and the remaining gluons have charge 0. 

Again, the total flux of a gluon inside the slab of thickness /d = ""^e 
on both sides of the loop as seen by the loop is ^g. The other half is lost 
on the loop. So the contribution of a fixed gluon species with non-zero 
charge is —1. As we suppose the gluons to be independent, the probability 
distribution for all the species inside the slab will factorize into single species 
distributions -P(/), / the number of gluons of that species inside the slab. 
Only the 2k{N — k) species with non-zero flux will contribute. Because of 
the factorization: 

/ \ 2k{N-k) 

{v,{L)) = ij2m{-iyj (4.9) 

and with the Poisson distribution for P{1) we get J2i^{^)(.~^y — exp— 2r, 
with / = n(T)lY)A{L) the average number of the gluon species in the slab. 
It follows from Eq. ()4.8() that the tension equals: 

Pk = 8lBn{T)k{N - k) . (4.10) 
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It is the dependence on the strength k of the loop, which is typical for 
the quasi-particle picture. We have checked in perturbation theory that 
deviations of this behaviour start to develop only to three loop order. Just 
above the transition we have no reason to trust the loop expansion. If 
the strong first order transition found at A'^ > 6 really implies a quasi- 
particle picture one should simulate the loop just above the transition for 
its A;-dependence. Paraphrasing ref. [26], such a behaviour just above Tc 
would suggest that the upper limit to the interface tension r^^rf will scale 
like A^, in accordance with a strong first order transition! To understand 
this argument, consider the complex plane with the possible phase of the 
Wilson line P = ex.pik^ on the unit circle. In the deconfined phase the 
effective potential has degenerate minima in the Z(N) vacua. At T = Tc 
one has another degenerate minimum in P = 0, the confined phase, which 
in coexistence with any of the Z(N) phases, will have the same tension r^^^ 
because of the obvious symmetry in the plane. A region of space with P = 1 
is separated from a region where P = expifc^ by a wall given by p^, for 
which the tension r^ N from ea. H4.10() . If just above the transition Tc the 
deconfined phase with P = starts to form in between the two Z(N) phases 
(so-called "wetting") we have 2rc,rf < pk- If wettting is to be true for all 
Z(N) interfaces the upper limit follows. If we know the Z(N) spin model, 
that falls in the same universality class as our gauge theory, we can check 
these statements quantitatively^. 



5. Magnetic screening mass and spatial Wilson loop 

Not only the force law between heavy electric charges like the heavy 
quark, but also the force between heavy magnetic charges tells us about 
the medium. The original idea of 't Hooft and Mandelstam was that of a 
dual superconductor, with the electric Cooper pairs replaced by some form 
of magnetic condensate. This condensate would be expected to screen the 
colour-magnetic field. 

In Section]^ we constructed an operator 14 (L) creating a magnetic fiux 
of strength expi/c^, Eq. (|4.5j) . This loop was space like. 

To get the monopole anti-monopole pair at points (0, r) we have the 
vortex end at and r on the positive z-axis. The vortex is given by a 
gauge transformation Vfc(x) which is discontinuous modulo a center group 
element exp ik^ when going around the vortex. The vortex is like the Dirac 
string in QED. It is unobservable by scattering with particles in the adjoint 
representation, as long as it has center group strength. It reads, up to a 
regular gauge transform: 



It is known that in the N state Potts model p^ — p — r^^d at the first order phase 
transition. 
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Vk = expiE.D{A)vkix,y) 



(5.1) 



with Vk{x, y) = arctan(|)y/j. When encirchng the point x = y = the gauge 
transformation expivkix,y) picks up a factor expz27rYfc = z^. This gauge 
transformation remains, by definition, unchanged along the z-direction and 
will be denoted by Vk{r). We say that 14 (r) creates a vortex or "Z(A'') 
Dirac string" of length r. That means, a Wilson loop W in the fundamental 
representation that encircles the vortex will pick up the Zk factor: 



VkWV,^ 



ZkW. 



(5.2) 



Any Wilson loop with non-zero A^-allity I will pick up a factor (zfc) • But 
Z(A^) neutral loops will not sense the Z(A^) Dirac string, hence the name. 

On the lattice the Hamiltonian operator will have magnetic plaquette 
operators. These will pick up the Zk factor and it is not hard to see that the 
Gibbs trace can be worked into a path integral along the usual lines, and 
on the lattice the latter takes the form: 



Fujr) ^ !DAexp-S(^k)iA) 
''''^ T ~ jDAexp-S{A) ' 



(5.3) 



The action St 



3(fc) is the usual lattice action, except for those plaquettes 
pierced by the Dirac string. Those plaquettes are multiplied by a factor 
expifc-^, as in Fig. El This string is repeated at every time slice between 
r = and r = 1/T. 



Fig. 6. Monopole antimonopole pair induced by twisting the plaquettes pierced by 
the Dirac string. 

Screening is expected in both confined and deconfined phases: 



Fuir) = Fuo -CM 



exp —mMr 



(5.4) 



All parameters are function of T. In the cold phase the screening is a 
consequence of the electric flux confinement. This is natural because the 
ground state contains a condensate of "magnetic Cooper pairs" , according to 
the dual superconductor analogy. It is a screening mechanism whose details 
are not understood. We dropped for notational reason the dependence on 
the strength k of the monopole in the coefficient cm- 
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In the hot phase there are indications from spatial Wilson loop sim- 
ulations that there is additional thermal screening from magnetic quasi- 
particles, as discussed in Section 121 The magnetic mass does not depend on 
the strength k of the source, just like the Debye mass does not, as discussed 
in the previous section. What we will show below is that the magnetic mass 
is a 0^"'' state in the fictitious Hamiltonian H. 

Analogous to the Wilson line correlator we consider the Hamiltonian 
H in the fictitious system of (x, y, r) variables. We search the operator 
Vk acting on the Hilbert space of physical states of this Hamiltonian, that 
reproduces the path integral Eq. 1)5. 3() ^ . So V^ should create a vortex in the 
(x, y) plane at every time slice r and the Hamiltonian H should propagate 
every one of these vortices in the z-direction over a distance r. So V^ is the 
't Hooft vortex operator discussed around Eq. (|5.1|1 : 

Vk = expi J TrE.D{A)vk{x,y). (5.5) 

with Vk{x,y) = arctan(|)Yfc. 

Both under parity (remember: only y -^ — y!) and charge conjugation 
the vortex Vk transforms into Vj. R parity leaves it invariant. Its spin J 
equals 0, despite the appearance of the rotated singularity line. On physical 
states the location of the singularity does not matter. Hence the operator 
ImVfc excites spin zero states with R = 1 and P = C = —1, and ReVfe 
excites scalars with R = P = C = 1. Which combination is the one we 
are after? In contrast to the Debye mass, we continue for T below Tc to 
measure the magnetic mass. At T = it becomes a three dimensional scalar. 
And with 3d rotational symmetry restored R and P parity are related by a 
rotation. Therefore, on a scalar state P = R. Thus ImV^ is excluded, and 
the magnetic mass is just that of the 0^"*"^. 

Perturbation theory is not reliable and we need lattice simulations [24]. 
Up to now these simulations are four dimensional and have limited accu- 
racy. They need to be repeated, also in dimensionally reduced form. Once 
they reproduce the mass levels of the fictitious Hamiltonian with sufficient 
accuracy, we can use them with more confidence for determining the tension 
of the space-like 't Hooft loops. 

At high T where reduced calculations are valid we expect to find the 
mass level of the 0"^"*" of the reduced Hamiltonian. From Teper's work [21] 



^ We use the same notation as for the vortex operator in {x, y, z) space as there is no 

risk for confusion. 
^ This fact was overlooked in earlier work, hep-ph/0308229 published in Vol. 40 of 

Erice International School of Physics, "From Quarks and Gluons to Quant umgravity, 

pg 72-73 and in Act. Phys.Pol. B, 2003, 12 : 5825. 
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the lowest 0^"^ mass in units of the string tension for SU(A^) {N > 3) gives 
in a large A^ expansion (see his table 34): 

^ =4.065(55)(l + 0.634(20)/A^2^O(l/iV^)). (5.6) 



For A^ = 2 one finds [21] for this ratio 4.718(43), for iV = 3 it is 4.329(41). 

5.1. The spatial Wilson loop 

The spatial Wilson loop is of interest because it monitors the magnetic 
activity in the plasma. At zero temperature it obeys an area law identical 
to that of time-like loop, controlled by the zero temperature string tension. 

The tension stays constant throughout the confined phase, and starts 
to rise about the critical temperature, indicating a new source of magnetic 
flux activity. 

Let us begin with some basics: a representation built from k fundamental 
representations is said to have A^-allity k. A center group transformation 
expz^ is mapped into expi/c^ in such a representation. Write a Wilson 
loop formed with such a representation as Wk (L) . Its average will then give 
an area law with tension ak- It is generally accepted that the Wilsonloop 
tension does not depend on the chosen representation, only on the N-allity 
k. 

At high T one can integrate out the hard modes, as they do not deter- 
mine the string tension. One can also integrate out the electrostatic modes, 
and wind up with a path integral controlled by the magnetostatic action: 

exp -ak{T)A{L) = — -^^ — -— — . 5.7 

The hard and electrostatic free energies /h and /e drop out in the ratio. 

The only dimensionful scale in the magnetostatic action is g^. So the 
tension, having dimension (mass) , can be written as: 

ak{T) = Ckgt,{l + Oig^)). (5.8) 

So the dominant contribution to the tension is entirely from the magne- 
tostatic sector. In Fig. you see a fit of the tension data to this parametric 
expression for SU(3). 

The authors took for the magnetic coupling g'^ = g'^, so neglected renor- 
malization effects of the scale gT, which are a few percent at T = 2Tc (see 
ea. ()3.7() ) On the other hand they included two loop renormalization effects 
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Fig. 7. The temperature over the square root of the spatial string tension CTs versus 
T/Tc for SU(3). The dashed Hne shows a fit according to a two loop scaling formula 
for the couphng, see text below Eq. H5.8|l . From Ref. [15]. 



in the coupling. Dropping the latter effects (as we argued before), and tak- 
ing into account the renormalization of qm as in ea. ()3.7() the data points in 
fig. constrain through a x^ analysis the values of Ck=i and Tc/Aj^o- We 
show this in fig.®. 

The value of Ck=i = 0.5530(20) follows from ref. [21] and is not compat- 
ible with the confidence limits. The central values for the ratios Tc/Aj^g 
are off the quoted values [28] . Taking only data points above 3Tc does not 
improve the situation. 

It follows that 4D data for higher than 4Tc are needed in order to see 
the 3D physics. For more detail the reader should consult ref. [28]. 

The spatial Wilson loop measures in a sense to be specified later the 
magnetic flux in the system. As we will see in section @ the 3D loops are 
well described by a magnetic quasi-particle picture. So this picture may be 
valid down to a few times the critical temperature. 

On the other hand the tension is flat from T = to T = Tc within errors, 
according to the data. In all of the confined phase the magnetic activity 
does not change. 
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Fig. 8. Confidence limits (90 resp. 68 % for outer, resp inner ellipse) for c^^i 
(abscis) and Tc/Ajj^ (ordinate) for SU(3) from the data down to 2Tc in fig.(I7|. 
FromRef. [281. 



6. A simple model 

In close analogy with the 't Hooft loop one can do a quasi-particle cal- 
culation for the k-tension of the Wilson loop. But what are the magnetic 
quasi-particles? 

We assume here that the magnetic screening length Im defined in section 
(0) defines three dimensional non-perturbative lumps, called magnetic quasi- 
particles. They are constituted by the magnetic gluons. 

Such a lump at rest will have a screened colour magnetic field. 

Note that the magnetic screening length defines a volume in which many 
elementary quanta, the magnetic gluons, are present, just like in the case of 
the Debye screening length. The difference is that the Debye sphere, though 
containing many elementary gluons , is not a bound state of the particles it 
contains. 

We are making the simplest possible assumptions: 

•The magnetic quasi-particles have a screening length Im ~ g^T much 
smaller than their average distance, i.e. they form a dilute gas. 
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•The magnetic quasi-particles are in the adjoint representation of the 
magnetic SU{N) group. 

The global magnetic group as introduced in ref. [19] should be either 
SU{N) (in case there is only Z(N) neutral matter) or SU{N)/Z(N) (in 
case Z(N) charged matter like quarks couples). So our adjoint representation 
covers both. 

We should have, like for the 't Hooft loop, a magnetic flux representation 
for the /c- Wilson loop. One can argue, starting from existing results [33], 
that for the irrep constructed by fully antisymmetrizing k quark representa- 
tions the average of the corresponding Wilson loop can be computed from: 

{Wk{L)) = (ew{igJdS-TrBYk}\ , (6.1) 

where, as before, Y^ = diag jj'iN — k,N — k,...,N — k,—k,... — k) with 

N — k entries —k and k entries N — k, so that it generates the center group 

element: 

27r 
exp (i27rYfc) = exp —ik — = z^. (6.2) 

For details on eq. (|6.1j) the reader should consult Appendic B and C. 

The Ifc-charge of a magnetic quasi-particle is it— with the same mul- 
tiplicity 2k(N — k) as for the electric quasi-particles. It contributes —1 to 
the Wilson loop Eq. (|6.1|) because only one-half of its flux goes through the 
loop. 

With our assumptions we can now, in precise analogy with the calcula- 
tion of the 't Hooft loop, understand why the /c-tension scales like: 

ak = ck{N -k)lMnM{l + 0{5)). (6.3) 

with 6 the diluteness parameter I^^um- This parameter, as argued at the 
end of the next subsection, turns out to be about 0.06. As in the gluon case, 
um is the density of one quasi-particle species. The constant c represents the 
effects of using the Poisson distribution for a single quasi-particle species and 
some elementary geometry, as in section 0, ea. (|4.10() ). So it is a number, 
in which colour dependence is absent. It is calculated in Appendix A to be 
c = 6.57128.... We will come back to this point in the next subsection. 

The simplest test of this k-scaling is to consider the ratios of tensions in 
which everything drops out, except the k-dependence. The tensions, as far 
as their leading term is concerned, can be simulated in 3 dimensions. 

Mathematically the degeneracy of the number of charged particles in 
the adjoint representation has the same k dependence as the quadratic 
Casimir operator of the totally antisymmetric representation built from 
k < N quarks. This quadratic Casimir equals CFk{N — k)/{N — 1), with 
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Cp = (N"^ — 1)/2N the the Casiniir of the fundamental representation. We 
will refer to our predicted k-dependence as Casimir scaling. 

6.1. Comparison to lattice results 

Our prediction is for the totally antisymmetric representation with N- 
allity k. 

At the time the model was proposed, there were only tensions known 
for SU(2) and SU(3) groups. There our prediction for the ratios is already 
implied by charge conjugation. 

Since the last three years lattice data on string tensions for N > 4 have 
been taken. The ratios found [22] for the totally antisymmetric irreps are 
close — within a percent for the central value — : 

SU(4) : CTa/o-i = 1.3548 ± 0.0064 Casimir : 1.3333 

SU(6) : CT2/CT1 = 1.6160 ± 0.0086 Casimir : 1.6000 

as/ai = 1.808 ±0.025 Casimir : 1.8000 

The results are that precise, that you see a two standard deviation, 
except for the second ratio of SU{6). This deviation is natural, since the 
diluteness of the magnetic quasi-particles is small, on the order of a couple 
of percent, as we will explain at the end of this subsection. So we expect 
corrections on that order to our ratios. 

There is a less precise determination of the ratio o^jcrx = 1.52 it 0.15 in 
SU{5) [23]. But the central value is within 1 to 2% of the predicted value 
3/2. 

The SU{8) ratios are known on a rather course lattice [23] and using a 
different algorithm: 



as/cTi = 1.692(29) Casimir : 1.714 
as/ai = 2.160(64) Casimir : 2.143 
a^/ai = 2.26(12) Casimir : 2.286 



In conclusion: the seven measured ratios are consistent with the quasi- 
particles being independent, as in a dilute gas. The number of quasi-particle 
species contributing to the k-tension is 2k{N — k). 

A posteriori, one can understand the diluteness being small. Take our 
formula for the tension, and multiply it with l\j to get the diluteness: 



liinu < liiCTk (6.4) 
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The diluteness has no parametric reason to be small. But from lattice 
data [21] one finds typically 0.05 '^ , as we discussed in section (jSJ. 



6.2. Partial pressure, density and N dependence of the tension 

There is an important question left. It concerns the behaviour of the 
tension for large N. In our expression for the tension of the Wilson loop 
we had a constant c, independent of the number of colours. The magnetic 
screening Im is 0(1) as argued above. Now, if we would take the density of 
one species to be: 

nM = d{N){cfNTf, (6.5) 

then the large N behaviour of cti would be determined by the product 
d{N)N . If we would take naively d{N) = 0(1), that is the density of 
one species to be 0(1), this would contradict the result from all orders in 
perturbation theory, that the Wilson loop tension is 0(1) for large A^! So 
we are forced to insist on a one species density of order ^, i.e. d{N) = -^. 

This dependence on N, combined with the idea that the magnetic quasi- 
particles are free, has obviously consequences for the pressure of the plasma 
discussed in sectional 

The effects from the magnetostatic action showed up at 0{g^). This 
contribution as discussed extensively there, contains the Stefan-Boltzmann 
factor ~ (iV^ — 1)T^, multiplied with {g'^N)'^{ and eventual non-leading 
terms of order ■^). 

If we believe - as seems reasonable from the lattice data discussed before 
- that the magnetic quasi-particles are independent, then they should deliver 
a partial pressure proportional to A^^ — 1 and the density of a given fixed 
species tim = jj'— — t^- We put the 't Hooft coupling in front of the 
temperature as it should in the 3d theory. 

Let the free energy per quasi-particle be f{N)T. So the partial free 
energy from all the quasi-particles will be (A^^ — ^)j^— — T^f{N)T, with 
the product ^f{N) of 0(1). 

So this implies f{N) ~ N for the free energy per magnetic quasi-particle. 
In other words the colour degrees of freedom of the lumps do show up in 
the free energy. It means the lumps have an energy A^e — 5, where 5 is the 
binding energy growing less fast than N . Any simulation of the partial free 
energy would then tell us the value of df{N). 



^ This small ratio is naturally explained in flux tube models of tlie spectrum. The 
electric analogue of this ratio is l^n ~ l/g^, and is large precisely for parametric 



Quasi- Particles in Hot QCD 31 



7. Conclusions 

In these lectures we concentrated on the quasi-particle picture of the 
electric and magnetic sector. The lattice data are consistent with the pre- 
dictions within a few percent, the typical order of magnitude of their dilute- 
ness. 

The magnetic quasi-particle model is certainly a viable idea, and it de- 
serves study in the Minkowski description. After all, this is where plasmas 
live. It is perhaps not superfluous to remind the reader that magnetic quasi- 
particles have a density that is suppressed by (log(;^))^. At asymptotic 
temperatures we have the usual Stefan-Boltzman gas. 

Of course, lattice simulations stay of invaluable support and are an in- 
dispensible testing ground. A most probably incomplete list of suggested 
simulations follows. 

• Determination of the partial three-dimensional pressure from lattice 
simulations [29]. It would further constrain the model. 

• At large N the strong first order results reported by Teper [26] suggest 
a quasi-particle picture down to the critical temperature and can be 
put into evidence by simulating at T^ spatial Wilson and 't Hooft 
loops. The simulation of the latter is now getting in a new stage [10], 
where we start to learn [27] the systematic errors from comparison 
with the known screening masses. The area law of the spatial 't Hooft 
loop permits simulation by multilevel methods [23]. From the results 
on the 't Hooft loop at T^ one could infer the corresponding Z{N) 
spin model. 

• If at large N indeed the quasi-particle picture is valid down to Tc, one 
would expect the topological susceptibility to drop to zero at Tc. This 
tendency is confirmed in a recent paper [30] and by the Pisa group [31]. 
Here one needs to distiguish exponential decay due to screening, valid 
for all N. 

On the theory side there are many questions that need more scrutiny. 

• Our model gives results for the AS irrep, depending on the N-ality. In- 
dependence of the representation is vindicated by lattice methods [32] 
for SU93). Taking this independence for granted there are argu- 
ments [34] that the large N limit of ak admits only corrections from 
powers of l/N"^. Our leading order result is inconsistent with this 
claim, but it may be that the series in the diluteness will undo the 
1/A^ corrections. Of foremost importance is the apparent accuracy of 
our result, better than 1/iV^ for A^ < 4. 
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• 



The use of the flux formula, as advocated in ref. [33], can be justified 
in a more specific way. In appendix C we argue why eq. (|6.1|1 is the 
one that applies to the AS irrep. 

• An amusing speculation can put the Wilson loop and the 't Hooft loop 
on an equal footing. Suppose that the magnetic group is operating 
at scales g^T as a local gauge group. Then we could define the Wil- 
son loop as a magnetic gauge transform with discontinuity in Z{N). 
The area law would then trivially depend only on N-allity. Also the 
magnetic flux formula for the loop, eq. 1)6. 1(1 . would be implied. 

• How do magnetic quasi-particles behave at high density, especially in 
the CFL phase with a VEV for the glue? Does one get bound state 
formation? 

It is amusing to make an estimate of the number of magnetic quasi- 
particles in the typical lattice volumes so far considered. This can be easily 
extracted from the values of Im in terms of the lattice length a (> 1), and 
the lattice size L in terms of a (~ 20). The diluteness of one species, eq. 
1)6. 4|) . gives a number on the order of a few hundred for one species and say 
iV<4. 

How can one view the mechanism of the transition in terms of the quasi- 
particles? First, at temperatures slightly above Tc, one is confronted with 
a thermal de Broglie wave length T of the same order as the interparticle 
distance g^T and the Bose-Einstein character of the quasi-particles starts 
coming into play. So the statistical independence is no longer a viable 
approximation. One expects the magnetic quasi-particles to condense in 
the groundstate below Tc- There they give rise to a non-zero string tension 
through electric flux tube formation. The simplest assumption-reasonable 
because also in 4D the diluteness as given by the ratio l\ja is small, about 
0.09- is that one has a condensate of our dilute Bose gas . To leading order 
in the density expansion, there are now however two and higher particle 
correlations. They will give deviations from the single particle multiplicity, 
that gave us in the plasma phase the Casimir scaling. 

At the same time their effect on the spatial Wilson loops is constant in 
T < Tc as we know from the simulations. This is an important clue to un- 
derstand. The tension of the time-like Wilson loops is related to the fraction 
of superfluid quasi-particles and is known to diminish as the temperature 
rises. So one would expect this tension to decrease when approaching Tc 
from below. This is qualitatively the case in lattice simulations. 

I thank Pierre Giovannangeli, Ben Grypaios, Alex Kovner, Laurent Lel- 
louch, Rob Pisarski, David Skinner, Mike Teper, Mikko Laine and Urs 
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Wenger for constructive and sometimes hard-needed deconstructive criti- 
cism. A Royal Society grant permitted me to write these notes in the 
stimulating atmosphere of the Oxford Dept of Theoretical Physics. Special 
thanks go to Harvey Meyer for communicating his preliminary SU{8) data 
on Wilson loops. 

8. Appendix A 

It is of interest to figure out how in our independent quasi-particle pic- 
ture the average of 14 (L) will behave as function of its strength k. As a 
byproduct we treat the geometry more precisely. 

A simple argument shows that it must be less than linear in k. First 
note that Vfc(L) = (Vi(-L)) . Imagine now we pull the k factors apart into 
an array of k parallel Vi (L) loops separated by distances well over twice the 
screening length. Then the loops get uncorrelated and the result for the 
array is k times the result for Vi{L). Thus we have a linear law, kpi, when 
the unit strength loops are far apart, because each of them gets disordered 
by the fluxes of the particles in its surrounding slab of thickness 2//) . 

Now we start to collapse the array. If we keep the loops in the array 
just one Debye length apart, any charge in a slab will disorder the two loops 
limiting the slab. This introduces correlations and the linearity may change. 

An estimate is based on the idea that infinitesimal slabs at distance I, 
and of width dl, are independently contributing to the average < Vk >■ So 
the average is built up from a product 11/ < Vk >i4h 

Let us denote by $(/) the flux normalized by the charge e at distance I 
from the loop, thought to be large and circular with radius R and area A. 

The average from one charge species from two slabs at distance it/ and 
of thickness dl is: 

< Vk >i,di< Vk >-i,di= Y. P{m)exp{ikm<^{l)} ^ P{n)exp{-ikn<^{l)}, 

m>0 n>0 

(8.1) 
where we used the independence. 

Take for P{n) the Poisson distribution to get: 

< Vk >i,di< Vk >^i,di= exp{-dlAn{T)2{l - cos{2kTT<^{l)))}. (8.2) 

So the contribution from all slabs is: 

< Vk >= exp{-2An(T) / dl{l - cos(2/c7r$))}. (8.3) 

Jo 

Note that a unit strength loop gets the same value when the gas is 
consisting of charges of strength k. 
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For the total flux going through our loop one finds in terms of the 
screened potential Aq = -^ exp —rjlD very simply: 



$(0 = — (Ao(r = 0-^o(r = V?2Ti?2)). (8.4) 

e 

Plugging this into eq. H8.3|) and dropping the corrections due to the size 
R of the loop, one gets for the tension for one single species: 

Pk{T)=2n{T) dl{l-cos{kTTexp-l/lD)). (8.5) 

Jo 

For k = 1 the integral /q°° dx{l — cos(7rexp —x) = 1.64282.... 
The derivative with k of p^ is easily computed: 

^p,(T) = 2^-^^{l-cos7Tk). (8.6) 

For odd k one finds from this: 

Pk{T)= 41 Dn{T) log k + p^{T). (8.7) 

The even k lie just below the logarithmic curve, and have derivative 
zero. So this k-dependence is like a logarithmic staircase. 

As we noted above, the Z-ionic tension for the unit fiux loop is identical 
to the electronic tension of a loop with Z units of fiux. So the slow rise of the 
tension with Z fiux units shows there is binding compared to a configuration 
of Z distant loops with unit fiux. 

Note that for the non-Abelian k-loops there is no k dependence for a 
single quasi-particle species, since the charge is ±(7, independent of k. The 
k dependence comes in through the multiplicity of the charged particles. 

9. Appendix B 

In this appendix we briefiy indicate the group theory needed to get 
from a given Young tableau (defining the irreducible representation R) the 
corresponding highest weight and the value of the quadratic Casimir.In what 
follows we suppose a representation to be irreducible without mentioning 
so. 

Let the Young tableau have ni boxes in the first row, 77-2 in the second 
row, etc.. Then one can define the non-negative numbers wi = ni — n/+i. 

Now the highest weight matrix for the Young tableau is defined from 
our Y matrices, defined above ea. ()6.2() by: 

N-l 

1=1 
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For example, for the totally antisymmetric tableau of k boxes in one 
column we have Hji = Y^. For the totally symmetric tableau with all k boxes 
in one row Hr = kLi. Note that the stability group of Y^ (the subgroup of 
SU{N) matrices commuting with Y^) is SU{k) x SU{N - k) x U{1). 

So the totally antisymmetric representation with k squares has a highest 
weight with this stability group. All other representations with k squares 
have different stability groups. 

We define one more diagonal NxN matrix by: 

N-l 

2y = 2 ^ = diagiN -1,N -3, ...., -N + 1). (9.2) 

1=1 

The quadratic Casimir operator C2{N,k,{wi} = C2{R) is defined by 
summing the square of all generators Ta in the representation R. The result 
is Z]^a = C2{R)Ir, where 1r is the unit matrix in R and C2{R) is a c- 
number (normahzation is [Ta,Th] = ifabcTcJabcfbcd = N5ad)- Then the 
quadratic Casimir equals: 

C2{R) = \{TrHl + 2TrYHR). (9.3) 

The quadratic Casimir for the fundamental representation is Cp = 
Ljy_j^N The Casimir for the antisymmetric representation is then 

C2{R = AS)=CF ^^^2iy (9-4) 

and for the symmetric representation it is: 

C2{R = SS) = c/-^±^. (9.5) 

One can show that for fixed k < N the antisymmetric Casimir is the 
minimal one. 



10. Appendic C 

In this appendix we review the Stokes formula [33] for the Wilson loop 
in the representation R. Let Vt be any gauge transformation that is periodic 
on the loop. Then, with V = 3 — igA: 

Wr{L) = f Dnexp{ig ( dS.Tr[HR{VLBVL^ - -VQ. x VJl^")]}. (10.1) 
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If R is the antisymmetric representation we have from the formulae for 
the highest weight in Appendix B the equation in the text, ea. H6.H) . Note 
that this representation is the unique one with stabihty group SU{k) x 
SUiN -k) xU{l). 

The integration over the regular gauges has been dropped in the main 
text because physical states do not feel them. Physical states include 
monopole configurations, which is why the second term in eq. (|1U.1I) is 
discarded in ea. H6.1|) . 

There is a comment related to this Stokes formula (see Appendix D). It is 
derived under the assumption that it is regulated by the SU(N) asymmetric 
top [33]). The question is whether the pure Yang- Mills theory average can 
be provided with such a regulator. For N = 2 and in three dimensions the 
answer is affirmative [35] by adding an adjoint Higgs system and letting 
the VEV go to zero, followed by decoupling the Higgs in the infinite mass 
limit. For general N the answer is analogous, but Nature realizes only a 
limited set of Higgs phases. They are limited to breakings of the type where 
SU{k) X SU{N -k) X U{1) is still unbroken, k integer and < [iV/2] [36]. 
That implies that the Stokes formula is only valid for those highest weights 
that have this symmetry, i.e. the totally antisymmetric ones with k boxes. 
This justifies eq. 1)6. 1|1 in the main text. 

11. Appendix D 

The derivation of ea. H10.1|) is based on Fourier analysis on the group 
SU{N), and on simple properties of the quantum-mechanical SU(N) rotator. 

We are interested in the Wilson line between two points x{si) and x{s2), 
the line between the two points being parametrized by s. The line is the 
ordered product of some irrep Rq with highest weight Hq of unitary matrices: 



W{s2,si) = Pexpig f ^ A,ds. (11.1) 

J Sl 



Here Ag = -^.A the projection of A on the line. 

W{s2,si) is a unitary matrix in the representation Rq with highest 
weight Hq. 

Along the line one can write the vector potential as a pure gauge 

A, = —UdsU^ (11.2) 

and so one can gauge away Ag except at the end points : 

W{s2,si) = D'^{U2Ul). (11.3) 
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To get the Wilsonloop we let S2 — > si and take the trace of eq. (|ll.:-{|) . 
In this form it is not hard to show that the Wilson loop is the Fourier 
component of the propagator of the SU(N) rotator with Hamiltonian 7i: 

f dn{nu2\exp{i{s2-si)n}\nui). (11.4) 

in the limit of vanishing moment of inertia /. To define the Hamiltonian 
we start with the Lagrangian L of the rotator. Define the angular velocities 
Va = TrQds^^Xa- In terms of the Cartan basis Xij,i ^ j = 1,....,N and 
Xd, d = 1, ..., A^ — 1 it is given by: 

^ = ^{E ^J + E Vd} + E VdHd. (11.5) 

ij d d 

The highest weight Hq of our representation Rq is written as Hq = 
J2d-^dXd- This Lagrangian gives the Hamiltonian in terms of the canonical 
momenta Ja- 

'H = jj{T.4+T.(J'i-H'i)'} = YAC2-J2J!+J2{jd-Hdf} (11.6) 

ij d d d 

Inserting a set of intermediate states, eigenstates of the Hamiltonian, 
labeled by -R,and {nid = J{R)d}, gives for the matrix element of the prop- 
agator {U2^2\ exp {—i{s2 — si)H}\UiQi): 

Y: {U2n\R,{md})exp{-z^-^^^{C2{R)-J2^d)} (11-7) 

R,{ma} d 

X exp{-i^^^jj^{J2{md-Hdf)}{R,{md}\Uin). 

d 

Now we have, writing m = m^: 

{U2^\R,m){RMUi^) = dRD^„^{U2mUi^)^) (11-8) 

with dji the dimensionality of the representation. Let / ^ in the expo- 
nent. This forces H^ = rrid and C2{R) = C2{Ro)- The other intermediate 
representations have higher quadratic Casimirs. 

So in this limit the integral over the matrix element, ea. H11.4|l . reduces 
to 

J dndn^D§',,{nU2{nUi)^)ew{i{s2 - si){C2{Ro) - ^TrHi)/2I}. 

(11.9) 
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Now the representation matrices D [Q) obey well known orthogonality re- 
lations: 

j dnDli{Q)D^^^{n^) = d~\R)6R,R,6k,n6i,m (11.10) 

Substitute eqn Ull.lUf) into the l.h.s. of ea. ()11.9|) to obtain the Wilson- 
loop as in ea. H11.3|) . multiplied with the phase factor 
exp{i(s2 - si)iC2iRo) - \TrHl)/2I]. 

The matrix element of the propagator can be written in path integral 
form 

/ Dfi(s)exp-i / dsL (H-H) 

JUiQ J 

with L as in ea. (|11.5j) . The final form of the Wilson loop in the irrep Rq 
and with weight Hq is then: 

W= f Di}exp{- i dsTrHo{n{V -igA).-^n^)} (11-12) 

The line integral on the r.h.s. can be easily transformed into a surface 
integral, which gives then eq. (|l().l)) . 
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